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X  ABSTRACT 

/  1  ■  ■  i 

An  analysis  of  solar  velocity  data  obtained  at  the 
Stanford  Solar  Observatory  shows  the  existence  of  solar  glo¬ 
bal  oscillations  in  the  range  45  to  105  microHz  (160  to  370 
minutes).  These  oscillations  are  interpreted  as  internal 
gravity  modes  of  degree  L* 1  and  ,1*2.  A  good  estimate  of  the 
order  of  the  modes  has  also  been  made. 
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Observations  of  the  global  solar  velocity  field  have 
been  recorded  at  the  Stanford  Solar  Observatory  since  1976 
(1).  Similar  recordings  from  the  Crimean  Astrophysical 
Observatory  are  available  starting  from  1974  (2).  A  spec¬ 
tral  analysis  of  the  combined  data  showing  the  prominent 
160.01  minute  peak  was  recently  published  (3).  A  careful 
examination  of  the  combined  Stanford-Crimea  spectrum  has 
convinced  us  that  there  may  be  a  few  other  peaks  in  the  low 
frequency  portion  of  the  spectrum  which  are  significantly 
above  the  noise  level.  For  that  reason,  we  decided  to 
reanalyze  the  original  Stanford  data. 

The  Data  Set  and  Analysis 

The  previous  report  described  the  observing  and  data 
reduction  procedures  in  detail (3).  The  observations  are 
differential  measures  of  the  line-of-sight  velocity  of  the 
solar  surface.  They  are  made  by  comparing  the  average 
doppler  shift  from  the  center  of  the  solar  disk  with  the 
doppler  shift  from  a  concentric  annulus.  Since  it  appeared 
that  interesting  features  were  present  in  the  low  frequency 
part  of  the  spectrum,  we  extended  the  computation  down  to 
v*o.  Since  only  the  low  frequency  part  of  the  spectrum  is 
examined  here,  the  data  was  averaged  into  5-minute  intervals 
and  normalized  within  each  day  as  in  the  previous  analysis. 
First  the  Fourier  transform  of  the  4  years  of  Stanford  data 
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(1977-1980)  was  computed  using  a  standard  fast  Fourier 
transform  code  (fft).  It  was  found  that  the  largest  power 
is  in  the  range  45  microHz  -  105  microHz  (about  160  to  370 
minutes) .  The  resulting  power  spectrum  shows  a  number  of 
sets  of  lines  with  separations  corresponding  to  day  side 
bands  (i.e.  11.57  microHz)  which  obviously  come  from  the 
nightly  gaps  in  the  data.  Comparing  the  4-year  spectrum 
with  individual  yearly  spectra,  we  discovered  that  the 
strongest  lines  were  more  prominent  in  1979  than  in  the 
other  years.  This  may  be  due  to  the  relatively  clear  skies 
and  the  distribution  of  observing  times  in  1979.  As  a  first 
step,  we  then  focus  our  analysis  on  that  year.  Figure  la 
shows  the  1979  spectrum.  The  observations  were  begun  on  7 
April  and  continued  through  23  July  with  most  of  the  data 
collected  in  late  Nay  through  July.  There  were  a  total  of 
240  hours  of  data  available.  Although  the  analysis  was  done 
using  daily  normalized  data,  the  results  were  rescaled  using 
the  average  normalizing  factor  to  provide  an  approximate 
scale  in  m/s. 

An  estimate  of  the  number  of  peaks  in  the  spectrum  that 
are  significantly  above  the  noise  can  be  made  by  examining 
the  cumulative  distribution  of  the  spectrum.  A  power  spec¬ 
trum  computed  from  a  normally  distributed  noise  source  will 
have  an  exponential  distribution,  thus  the  logarithm  of  the 
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cumulative  distribution  will  decrease  linearly  with  power 
(e.g.(4)).  A  departure  from  a  straight  line  is  an  indica¬ 
tion  of  the  presence  of  significant  spectral  features  and 

2 

the  slope  of  the  line  is  a  measure  of  the  variance  (or  )  in 
the  spectrum. 

Figure  2  shows  this  plot  for  the  45-105  microHz  range 
of  the  1979  spectrum.  The  actual  data  is  represented  by 
large  dots.  The  smaller  points  are  for  a  spectrum  computed 
the  same  way  from  data  constructed  by  taking  the  original 
data  for  each  day  in  reversed  order,  thus  keeping  the  origi¬ 
nal  window  function.  Any  coherent  signal  present  in  the 
original  data  will  be  eliminated  by  this  procedure.  The 
noise  statistics  for  the  modified  data  power  spectrum  should 
be  unchanged  for  periods  up  to  the  average  daily  observing 
time,  i.e.  for  frequencies  higher  than  40microHz.  It  can  be 

seen  that  there  are  more  peaks  with  value  above  0.20 

2 

(m/s)  /freq  than  in  the  reversed  data.  The  level  at  which 
the  actual  distribution  begins  to  depart  from  the  noise  dis¬ 
tribution  is  around  2.5  or.  Since  the  spectrum  was  computed 
with  a  resolution  of  0.02  microHz  but  has  a  natural  resolu¬ 
tion  of  only  0.11  microHz,  there  will  be  about  five  points 
shown  for  each  significant  peak  in  the  spectrum.  Also, 
since  the  data  has  gaps  at  night,  the  day  sidelobe  structure 
introduces  2  to  4  apparently  significant  artificial  peaks 
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for  each  true  peak.  These  considerations  suggest  that  the 
largest  10  or  so  peaks  are  most  likely  not  noise. 

Since  the  spectrum  appears  to  be  dominated  by  the 
sidelobe  structure  from  the  observing  times,  a  procedure 
must  be  performed  to  eliminate  these  sidelobes.  An  itera¬ 
tive  peak  removal  technique  was  used  to  find  and  remove  the 
peaks  one  at  a  time.  First  an  fft  was  computed  and  the 
largest  peak  determined.  Next  that  peak  was  accurately 
found  with  a  fine  resolution  simple  Fourier  transform  in  the 
vicinity  of  the  peak.  Finally  the  corresponding  sinusoidal 
signal  was  subtracted  from  the  original  data.  This  pro¬ 
cedure  was  repeated,  producing  a  list  of  frequencies  free  of 
day  sidelobes.  The  list  of  the  first  22  peaks  found  in  this 
way  includes  14  peaks  in  the  45-105  microHz  range.  The  list 
of  14  peaks  is  given  in  Table  1. 

To  see  the  effect  of  removing  these  peaks  from  the  ori¬ 
ginal  data,  the  spectrum  of  the  residual  data  after  sub¬ 
tracting  the  22  sinusoids  is  shown  in  figure  lb.  Figure  lb 
is  shown  to  the  same  scale  as  figure  la.  It  is  clear  that 
far  more  than  14  peaks  have  been  removed  in  the  resulting 
spectrum.  This  method  of  peak  identification  appears  to  be 
a  powerful  tool  for  analysis  of  power  spectra  computed  from 
data  with  complicated  windows.  The  method  finds  a  table  of 
periods,  phases  and  amplitudes  that  best  clean  the  spectrum. 


Table  1 

Frequencies  of  17  peaks  found  in  the  45-105  microHz  range  of 
the  1979  spectrum.  The  frequencies  are  in  microHz  with  the 
corresponding  periods  shown  in  minutes.  The  power  is  shown 
as  (m/s)  per  0.11  microHz.  The  classification  of  n  and  l 
are  described  in  the  text.  The  mode  in  parentheses  is  an 
alternate  identification  for  the  peak  at  62.30microHz  and 
its  associated  side- lobes.  Note  that  the  peak  at 
59.52microHz  could  be  identified  as  part  of  an  1*1  or  1*2 
series.  The  modes  marked  with  an  asterisk  are  also  found  in 
the  spectrum  from  the  Crimean  Astrophysical  Observatory (9) . 
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Prom  experiments  with  artificial  data,  we  have  concluded 
that  while  the  correct  peak  is  usually  found,  a  side  lobe 
introduced  by  the  window  is  occasionally  selected.  This  is 
a  weakness  of  the  procedure  that  requires  care  in  examining 
the  resulting  peaks.  We  have  also  found  that  the  amplitudes 
determined  from  data  with  a  complicated  window  are  not  to  be 
relied  upon. 

The  method  described  above  relies  upon  the  single 
assumption  that  the  examined  signals  consist  of  coherent 
sine  waves.  It  is  then  not  surprising  that  for  each  fre¬ 
quency  removed,  a  complete  day-sidelobes  set  of  lines  disap¬ 
pears  from  the  spectrum.  A  further  check  is  necessary  to 
determine  whether  or  not  the  signal  is  coherent  for  the  full 
span  of  the  observations.  For  that  purpose,  the  data  with 
the  tabulated  frequencies  removed  was  divided  into  two 
parts,  7  April  through  15  June  and  16  June  through  23  July. 
The  spectra  of  each  part  separately  did  not  show  any  of  the 
removed  frequencies.  If  any  of  the  frequencies  removed  were 
not  present  in  both  halves  of  the  data,  they  would  show  in 
the  spectrum  of  one  or  both  halves  of  the  cleaned  data. 

Two  of  the  fourteen  peaks  are  apparently  day  aliases  at 
1/4  and  1/8  of  a  day.  (With  only  3  months  of  data  the  pre¬ 
viously  reported  peak  at  160.01  minutes  can  not  be  dis¬ 
tinguished  from  1/9  of  a  day  and  is  masked  by  the  sidelobes 
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of  the  lrwer  day-harmonics) .  This  leaves  12  peaks  for 
further  analysis.  There  is  some  ambiguity  in  the  true  iden¬ 
tification  of  the  largest  of  the  remaining  peaks.  The  first 
1/day  sidelobe  is  of  the  same  amplitude  as  the  main  peak. 
This  could  be  due  to  two  true  peaks  with  separation  1/day  or 
2/day.  By  careful  examination  of  the  peak  shapes,  we  have 
identified  the  peak  as  a  large  peak  at  62.3microHz  with  a 
smaller  peak  at  96.9microHz,  although  the  identification 
could  have  been  made  as  a  true  peak  at  73.8microHz  with 
smaller  peaks  at  50.7  and  96.9microHz.  The  alternate  iden¬ 
tifications  are  shown  in  parentheses  in  Table  1. 


Interpretation 


In  order  to  interpret  the  peaks  found  in  this  spectrum, 
we  must  seek  guidance  from  theory.  Estimates  for  the  spac¬ 
ing  of  g-mode  oscillations  made  from  the  full  asymptotic 
approximation  (5)  and  calculated  from  complete  solar  models 
(6)  suggest  that  the  g-mode  oscillations  of  the  same  degree 
l  should  be  about  equally  spaced  in  period  for  order  large 
enough  (n  >  6).  Therefore  we  checked  the  list  of  prominent 
peaks  for  equal  spacing  in  period.  Three  of  the  four  larg¬ 
est  peaks  are  separated  in  period  by  15.5  minutes.  This  is 
very  close  to  the  spacing  of  the  standard  solar  model  for 
g-modes  with  degree  two  and  order  greater  than  10. 
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To  aid  in  further  identification  of  the  inodes  we  can 
use  the  asymptotic  approximation  (5)  which  shows  that  for 
sufficiently  large  order  n  the  period  T  is: 

T  =  To  (n  +  f  -  j)  /  \] Ul+  1 ) 

If  we  assume  that  the  three  peaks  with  15.5  minute  separa¬ 
tion  are  part  of  an  £=2  series,  we  can  calculate  T  and  then 

o 

the  probable  spacing  for  the  £=1  and  £  = 3  series.  Doing  this 
we  see  that  most  of  the  remaining  peaks  are  likely  to  be 
part  of  the  £-1  series.  If  we  have  an  estimate  for  £,  we 
can  plot  the  observed  periods  as 

T  \J £( £+1)  vs  (k  +■  1/2  -  1/4)  where  k  is  an  integer  increas¬ 
ing  with  period.  If  the  peaks  are  consistent  with  the  model 
and  we  have  assigned  the  correct  degree  £,  the  points  will 
all  lie  on  a  straight  line.  From  the  intercept  of  this  line 
we  can  determine  both  the  best  value  for  n  and  check  the 
-1/4  term.  This  has  been  done  in  Figure  3. 

The  three  largest  peaks  were  used  to  define  the  1=2 

series  and  thus  to  find  T  .  The  filled  circles  in  Figure  3 

o 

represent  the  three  largest  peaks  which  determine  the  1=2 
series.  The  small  dots  represent  all  the  other  modes  iden¬ 
tified  in  table  1.  The  line  was  found  from  the  three  large 
1=2  peaks  only.  It  can  be  seen  that  the  other  peaks  are 
organized  in  the  expected  g-mode  structure.  Mote  that  in 
several  cases  two  peaks  have  been  assigned  to  one  order  n  in 
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the  l- 1  series.  These  peaks  have  an  average  separation  of 
1.22microHz  and  could  be  evidence  of  rotationally  split 
modes.  The  Tq  implied  by  the  identified  modes  is  38.6+0.5 
minutes.  The  intercept  of  the  best  fit  line  with  the 
abscissa  is  expected  to  be  0.25  and  is  found  to  be  0. 2+0.1 
thus  the  order  n  is  most  likely  correctly  determined. 

Now  we  must  ask  if  the  Stanford  instrument  is  sensitive 
to  degree  l  -2  g-raodes  oscillations.  For  the  previously 
reported  acoustic  mode  oscillations  in  the  5-minute  range, 
the  Stanford  velocity  differencing  scheme  is  most  sensitive 
to  modes  of  degree  For  low  frequency  g-modes,  Gough  and 
Cristensen-Dalsgaard  (8)  have  shown  that  for  periods  around 
160  minutes  the  Stanford  instrument  is  most  sensitive  to 
modes  with  degree  l**5  to  £=7.  For  longer  periods  however, 
the  instrument  begins  to  be  relatively  more  sensitive  to 
modes  with  lower  degree.  This  variation  with  period,  unlike 
the  case  for  the  p-modes,  is  due  to  the  importance  of  the 
horizontal  motions  induced  by  the  internal  gravity  modes. 

The  largest  peak  at  62.29microHz  appears  not  to  be  con¬ 
sistent  with  the  asymptotic  formula  with  the  Tq  we  have 
found.  If  the  alternate  choice  of  the  62. 20microHz- 
73.84microHz  pair  were  chosen,  the  peak  at  73.84microHz 
which  would  be  identified  as£**l,  n«8.  In  this  case  a  peak 
at  50.7microHz  shows  up  in  the  peak  finding  procedure  and 
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would  be  identified  as  4»2,  n»20.  We  also  note  that  the 
frequency  of  this  peak  is  almost  identical  with  the  fre¬ 
quency  difference  between  5-minute  modes  with  order  and 
degree  (n,f»l),  (n,£»0)  and  close  to  the  difference  between 
(n,£»2),  (n,f*l).  This  beat  frequency  could  show  up  in  the 
Stanford  observations  through  a  non-linear  coupling  in  the 
Sun.  We  do  not,  however,  find  any  peak  near  136microHz 
which  would  correspond  to  beats  between  acoustic  modes  of 
the  same  degree  but  with  order  differing  by  one. 


The  relative  peak  sizes  and  simplicity  of  the  analysis 
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